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Abstract 

In the kinematic region of small x and large in deep inelastic scattering, 
presently being explored by HERA, we present an analysis of the evolution of the 
longitudinal structure function Fj^{x, Q^) in the double scaling limit of Ball and Forte. 
We fit the evolution to a 1/x^ behaviour and extract the value of A. We also study 
the behaviour oi R = F2/2XF1 — 1. We present comparisons of both Fl and R with 
the corresponding MRS fits in this region of x and Q^. 
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The recent spurt of activity in the physics of low x QCD has been guided by the electron- 
proton scattering measurements at HERA which have studied the proton structure function 
F2{x, Q'^) in the very low x region, x < 10~^. The steep rise of F2 at these small values of 
X have given rise to a flurry of theoretical activity in an attempt to understand the nature 
of this rise through 'standard' QCD. 

There have essentially been two approaches to the study of the HERA data. In order 
to study the effects of summing a^ln^ unaccompanied by InQ^, the preferred approach 
has been to study the BFKL equation which generates a singular behaviour for the 
unintegrated gluon distribution /(x, fcf.), with A = as41n2 (for fixed, not running a^). For 
running a^, A ~ 0.5. [Q. From this an asymptotic form for Fi{i = 2, L) can be inferred. 

The other distinct approach is to attempt to describe the data through a 'standard' 
Altarelli-Parisi (or DGLAP) evolution equation to the next-to-leading order approximation. 
Here again, the data imply a steep gluon distribution with the gluon density rising sharply 
as X decreases, even for comparatively low values of Q^. Ball and Forte [Q], in a series of 
papers, have used this approach to exhibit the scahng properties at low x, generated by 
QCD effects and have shown that the HERA measurement of F2{x,Q'^) is well explained 
by this approach. 

We describe their result very briefly. They show that at sufficiently small x and large 
t = In jj, the structure function F2 exhibits double scaling in the variables 



a 




where the starting scale for Ql in to = \n-j^ can be just a little more than Ql = IGeV'^, 
and A = Aqcd- Double asymptotic scaling results from the use of the operater product ex- 
pansion and the renormalisation group at leading (and next-to-leading) order, and predicts 
the rise of F2 on the basis of purely perturbative QCD evolution. 

The asymptotic behaviour in this approach (to leading order) has the form 

F^{a,p) ~ iV/(2)2^exp[27p-5(-)] 



where 

7 

and 



x[l + 0(i)] (2) 




2nf , ,^ 2nt 



^=(1 + TT1^)(1- 



with / an unknown distribution that depends on the starting value of the distribution. 

In this paper, we concentrate on the Ball and Forte approach to QCD evolution to 
study the other independent proton structure function - the longitudinal structure function 
for the proton, Fl{x, Q"^). F^ has not yet been measured at HERA and it will of course be 
very interesting to see how the data compares with theory. 
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Before we embark on the calculation, we would like to make a few comments about these 
different approaches. As is well known, in the traditional DGLAP or Ball-Forte analysis, 
there are two unknowns viz. F2{x,Q'^) and gi^x^Q"^) (the gluon distribution) and these are 
determined by two measurements viz. a direct measurement of F2 (from the cross section) 
and the variation of F2 which gives g{x,Q'^). These two measurements fix these two 
unknown quantities exactly since there are no free parameters in this approach that need 
to be fitted. In contrast, in the BFKL approach, due to the resummation of the In ^ leading 
to evolution also in x (by some exponent), there is a third parameter which needs to be 
fitted. Since there have been only two measurements till now, this parameter, in the BFKL 
approach, can only be fixed with a third measuremnent, for example, that of Fl. This is in 
striking contrast to the Ball-Forte DGLAP approach where there are no parameters that 
can be tuned in order to fit any future measurement of Fl{x,Q'^). A theoretical analysis 
of the behaviour of Fl and subsequent comparison with experiment is, therefore, the need 
of the hour. We now present our analysis of Fl{x, Q^). 

In perturbative QCD, the nucleon structure functions are determined in terms of the 
quark and gluon distribution functions which in turn follow the DGLAP equation, 

f.|p«G(Q^) (3) 

where, 

P= (p'' and G(x,g2) = 

are the splitting function and the parton distribution matrices respectively. The nonsinglet 
quark distribution follows a similar equation but we shall ignore it as the non singlet 
contribution to the longitudinal structure function is negligible. 

The above equation (3) can be solved by using standard Mellin transformation tech- 
niques [0] and in the asymptotic limit cr — > 00 yields the result (2) for the structure function 
F2. 

Note however, that in the standard case, the quark distribution functions are defined 
in terms of the structure functions F2, viz., 

F2(x,Q2) = M ^enx(g,(x,t) + qNs{x,t)) . (4) 

We could, however, start by using Fi as the defining structure function instead of F2 and 
then the definition of the former would yield 




which can be expressed in terms of the singlet and non singlet distributions as in the 
usual case, and where the superscript (1) denotes the fact that the distributions have 
been defined with respect to the structure function Fi. If we take the quark and gluon 
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distributions defined with respect to F2 as the standard distribution, we can compute q^^^ 
and q^^^ with respect to these and the resuh is 

g«(x,g2) = l^'^{q(y,Q^)[S(l - z) + a.A/f(;.)] + g{y,Q')a,Affiz)} (6) 

and a similar equation for q^^\ where we have 

z = - ; AfHz) = ; A/f = - ^2-^ 

y Sir 7T 

Now if we calculate the structure function Fi, in terms of these quark and gluon densities, 
using (5), the result, in terms of the standard parton distributions obtained from F2, is 
given by 0] 




=<?.<? 



HI - z) + ^al^{z) 



+ 9{y) 




27r 

(7) 

where and g are as usual the quark and gluon parton densities, Cj is the charge of qi and 
in analogy with F2, we have written the expression in terms of and aj^^ which are the 
transverse virtual photon-quark (antiquark) and the virtual photon - gluon cross sections 
respectively. The longitudinal structure function Fl = F2 — 2xFi can now be easily 
evaluated with the standard expression for F2, in terms of the quark and gluon densities, 
and, rewritten in terms of F2, yields the Altarelli-Martinelli formula, 

Fl{x,Q ) = X —i-F2{y,Q ) + —yg[y,Q ) I 1 

In deriving the above formula, one uses the lowest order result for Fi, by introducing the 
running coupling constant, as(Q^), and using the momentum dependent parton densities. 
This formula can be used to compute the value of Fl for low x and large Q^, which combined 
with calculated values of F2 will give an estimate of R. 

In the present work, we work in the "double scaling" regime, i.e we consider the asymp- 
totic solutions of F2 and g in the limit cr ^ 00 and fixed but large values of p where a and 
p are defined in (1). 

The evaluation of F2 and the gluon distribution function g and the fit of F2 to experi- 
mental values have been carried out in and we have used the same best fit normailiza- 
tions as in p|. In the absence of any existing data for F^, in this kenematical regime, we 
have fixed an overall normalization for Fl as 0.14 in order to closely match the various 
parametrizations for the same. Indeed, it would be possible to carry out a detailed anal- 
ysis of the possible boundary conditions along the lines of [0, once experimental values 
for Fl are obtained. We must mention here that we are interested in only the asymptotic 
form of the evolution of the various functions and the DGLAP evolution equations cannot 
fix the numbers absolutely and these must be uniquely determined by the experimental 
data. In particular, in a recent paper, De Roeck et al. have fitted the expression (2) 
for F2 to the latest measurement of the proton structure function by the HI collaboration 
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at HERA. We have assumed the result of their fit, which gives A = Aqcd = 248MeV and 
Ql = 1.12 GeV\ 

We have numerically evaluated equations (7) and (8) and subsequently Fl and R, in 
the double scaling limit. Since HERA has not measured either Fl or R, we have given a 
comparison to the MRSA-fit 0. The data for Fl is shown in figure (la) and (lb), while 
that for R is shown in fig (2). The fits are shown on the same graph. As is clear, apart 
from overall normalization, about the arbitrariness of which we have already commented 
before, the behaviour of Fl and R closely follows the fit. 

In order to get a somewhat more quantitative handle on the data, we have done a fit 
of the graph for Fl to the form a{Q'^)x~^^^ \ The result of the fit is shown in the Table. 
It is strikingly clear from the table that a((5^) reaches an asymptotic value of 0.0034 and, 
more interestingly, A(Q^) reaches a asymptotic value of 0.55 We will have more to say on 
this later. 

A few words about our results are now in order. In this paper, we have considered the 
form of the longitudinal structure function at very low values of x (10~^ < x < 10~^) for a 
wide range of Q^. We have also calculated R and compared it with the present MRSA-fit. 

We see clearly that the steep rise of the longitudinal structure function as predicted by 
hard QCD processes gives a behaviour of the form This behaviour is similar to that 

predicted by the BFKL [|l| analysis. In the BFKL analysis, summing the aslog{l/x) terms 
unaccompanied by logQ^ generates a singular x~'^ behaviour for the unintegrated gluon 
distribution. For a running coupling, the BFKL equation has to be solved numerically 
under some reasonable assumptions, and it yields A ~ 0.5. (For a somewhat detailed 
discussion, see, for example, |^). This behaviour subsequently gives a similar behaviour 
for F2 and Fl- 

The rise of Fl is very close to that predicted by the MRSA-fit even though the overall 
normalization is presumably different and will be a function of the input distribution. 

At this level of analysis, it is impossible to distinguish, purely from measurements of 
Fli the behaviour as predicted by BFKL and that of the Ball- Forte analysis. Thus the 
role of the hard pomeron in deciding the nature of the rise of the structure function is not 
at all clear, at this level of analysis. It is possible that very precise meaurements of Fl 
coupled with the non-leading contribution to the Ball-Forte result will allow us to make a 
clear statement on the role of the hard pomeron. However, if indeed, Fl is seen to go as 
x~^'^ asymptotically, while experimentally F2 goes as a;""'^, then, even though the absolute 
value of F2 at these values of x and is much larger that Fl-, it is possible to find some 
X value at which Fl > F2, if indeed this behaviour persists. If that happens then since 
2xFi = F2 — Fl, this would imply a negative LHS. Since the LHS is the momentum density 
and hence positive, some new physics must take over at this stage to reduce the rise of Fl 
to keep the RHS positive. 

A measurement of the longitudinal structure function and R will shed light on the 
importance of hard processes in the kinematic region being explored at HERA and will 
consequently give us more information about the gluon. 
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Table 









20 


0.0183 


-0.316 


50 


0.00777 


-0.443 


100 


0.00652 


-0.467 


400 


0.00483 


-0.505 


800 


0.00424 


-0.521 


1200 


0.00395 


-0.529 


1600 


0.00376 


-0.535 


2000 


0.00362 


-0.539 


2400 


0.00351 


-0.542 


2800 


0.00342 


-0.545 
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Figure Captions 



la and lb: Plot of Fl vs. x for different Q^. The dashed hnes are the MRS-fit shown for 

comparison. 

2: Plot oi Rys X for different Q^. The dashed lines are the MRS-fits shown for comparison. 



Table Caption 



Values of a((5^) and A((5^) as a function of Q"^ (in the fit for Fl ~ a{Q'^)x^). 
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Figure 1a 




Figure 1b 




Figure 2 



